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CAN THE LYMAN CONTINUUM LEAKED OUT OF H II REGIONS EXPLAIN DIFFUSE IONIZED GAS? 

KWANG-lL SEON^ 

Abstract 

We present an attempt to explain the diffuse Ha emission of a face-on galaxy M 51 with the "standard" 
photoionization model, in which the Lyman continuum (Lyc) escaping from H II regions propagates large 
distances into the diffuse interstellar medium (ISM). The diffuse Ha emission of M 51 is analyzed using thin 
slab models and exponential disk models in the context of the "on-the-spot" approximation. The scale height 
of the ionized gas needed to explain the diffuse Ha emission with the scenario is found to be of the order of 
~ 1 - 2 kpc, consistent with those of our Galaxy and edge-on galaxies. The model also provides a vertical 
profile, when the galaxy is viewed edge-on, consisting of two-exponential components. However, it is found 
that an incredibly low absorption coefficient of k,o ^ 0.4-0.8 kpc~^ at the galactic plane, or, equivalently, an 
effective cross-section as low as aeff ~ 10~^ of the photoionization cross-section at 912A is required to allow 
the stellar Lyc photons to travel through the H I disk. Such a low absorption coefficient is out of accord with 
the properties of the ISM. Furthermore, we found that even the model that has the DIG phase only and no H I 
gas phase shows highly concentrated Ha emissions around H II regions, and can account for only < 26% of 
the Ha luminosity of the DIG. This result places a strong constraint on the ionizing source of the DIG. We also 
report that the Ha intensity distribution functions not only of the DIG, but also of H II regions in M 5 1 , appear 
to be lognormal. 

Subject headings: galaxies: individual: M 51 — galaxies: ISM — ISM: H II regions — ISM: structure 



1. INTRODUCTION 

Diffuse ionized gas (DIG) is a major component of 
the interstellar med ium ( ISM) in galaxies (Revnolds 1991; 

; IWalterbos & BraunI Il996l: iHaffner et al.l l2QQ9h ! The large 
energy requirement of the DIG strongly s uggests that OB 
stars are the only viable ionization source (lReynoldslll984l) . 
The "standard" photoionization model of the DIG assumes 
that the ionizing photons that leaked out of H II regions 

■ in the galactic disk ionize the diffu s e ISM (Mathis 1986; 

I Sokolowski & Bland-Hawtho rnI 119911; iDomgorgen & Mathis' 

. 11994)) . If the Lyman continuum (Lyc, ionizing continuum) 
photons escaping from the traditional H II regions were the 
dominant source responsible for ionizing the DIG, then the 
Ha intensity from the DIG would be ex pected to co rrelate 
with that of discrete H II regions. Ferguso n et al.l (Il996|) found 
a correlation between the DIG and bright H II regions over 
both small and large scales. 

However, it is not clear how a major fraction of Lyc 
photons can escape from the immediate surroundings of 
an OB association and how the escaped photons can 

I be transmitted across di stances of the order of a kpc. 
iNorman & Ikeuchil (|1989|) suggested that superbubbles or 
"chimneys" around OB associations could i onize the halo 
gas above t he chimn eys. iMiller & Cox' (ll 993h . lDove & Shulll 
Jl994), and lDove. Shull. & Ferrara (2000) quantitatively in- 
vestigated the possibility of penetration of the Lyc into 
the Galactic halo. These senarios may explain the dif- 
fuse Ha emission, in particular vertical Ha filaments or 
"wor ms," seen in external edge-on galaxies, such as NGC 
891 (iRand. Kulkarni & Hesteilll990l: lDettm^ll990h . How- 
ever, the Lyc must travel not only large distances above the 
plane, but also within a galactic plane (Reynolds 1993). This 
requirement seems evident when the diffuse Ha emission far 
from the classical H II regions in face-on galaxies are exam- 
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ined. 

IZurita et aH (120021) were the first to attempt modeling the 
global morphology of the DIG of a face-on galaxy with the 
Lyc escaping from bright H II regions. They assumed a thin 
slab disk with a constant absorption coefficient over a galaxy 
NGC 157 and compared the predicted morphologies with the 
observed data by varying the absorption coefficient. In this 
way, they found that the diffuse Ha surface brightness distri- 
bution is, surprisingly, well reproduced by assuming no sig- 
nificant absorption of ionizing photons. 

However, they neither compared the absolute Ha flux of the 
models with the observed data nor discussed the implications 
of the absorption coefficient obtained to explain the observed 
morphology. If OB stars are to provide the diffuse ioniza- 
tion, then the ISM must be not only transparent to ionizing 
photons, but also sufficiently dense to provide the observed 
Ha emission (Miller & Cox 1993). In the extreme case of no 
absorption, no Ha radiation is emitted. Therefore, it is nec- 
essary to compare the absolute Ha flux as well as the mor- 
phology of the Ha image. It is also worthwhile to consider 
a geometrically thick disk in the modeling and compare the 
scale height obtained to account for the observed Ha lumi- 
nosity with those of edge-on galaxies. 

The transmission of the Lyc into the DIG has been at- 
tributed to the inhomogeneity of the ISM, such as a frac- 
tal H I distribution produced by ISM turbulence ('E lmegreenI 
119971: IWood & MathisI 120041) . However, no attempt has been 
made to quantify the ISM density structure required to al- 
low the transmission of the Lyc. It is not clear whether such 
inhomogeneity of the ISM in typical spiral galaxies is con- 
sisitent with the requirement of the "photoionization" sce- 
nario. Meanwhile, it is known that the probability distri- 
bution functions (PDFs) of not only the local densities, but 
also the column densities of the turbulent ISM, are lognormal 
dYazquez-Semadeni 1994';'Qstriker, Stone, & Gammie"2005 
iFischera & Dopi ta 2004). With this property, it is possible to 
quantitatively access how turbulent or porous the ISM should 
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be for the stellar Lyc radiation to be responsible for the dif- 
fuse Ha emission with the "standard" photoionization model 
or whether the absorption coefficient obtained in the modeling 
is in accord with the current knowledge of the ISM. 

In this paper, we calculate thin slab models similar to those 
of Zurita et al. ( 2002) and more elaborated models incorpora- 
tiong an exponential disk for a face-on galaxy M 5 1 . We find 
that the "standard" scenario requires absorption too unrealis- 
tically small to be believed, but the obtained scale-height of 
the galactic disk is consistent with those of edge-on galaxies. 
We also report that the PDFs of the Ha intensities of the DIG 
and H II regions in the galaxy M 51 are lognormal. 

2. Ha DATA 

For our purposes, a face-on galaxy is the best choise since 
we can easily identify H II regions, i.e., the ionizing source 
that is expected to be the dominant source. When an Ha 
emission map is ready, the first step for the analysis is then 
to separate H II regions from the diffuse Ha emission. To 
do so, we used an automated photometry procedure known as 
"Hllphot," described in Thilker, Braun, & Walterbos (2000), 
in which they applied the procedure to the Ha data of M 5 1 . 
Therefore, as a matter of convenience, we use the same Ha 
data of M 5 1 as they did. Their publicly available data, pro- 
vided together with the "Hllphot" procedure, has the same 
angular resolution, but a sligh tly smaller field size than that of 
iThilker. Braun . & Walterbos (2000). 

For th e present analysis, we assurn ed a distance of 9.6 Mpc 
to M 5 1 (ISandage & TammannlI974l) . At this distance, a pixel 
in the image corresponds to a linear resolution of Ar = 32 
pc. No correction was attempted for extinction to compare 
with earlier studies. We detected 1356 H II regions. Of the 
total sample, 1061 regions were classified as "photometric- 
quality" detections having S/N > 5. The extent of each H II re- 
gion was determined usi ng a terminal surface brightne s s slop e 
of 1.5 EM pc-\ as in IThilker. Braun. & Walterb"osl (l2000h . 
Here, EM is an emission measure (cm~^ pc). A continuum- 
subtracted Ha image, where H II regions were masked, and 
the Ha intensity histograms for the DIG and H II regions are 
shown in Figure [T] It is evident in the figure that the Ha in- 
tensity histograms are lognormal, both in the DIG and H II 
regions. The width of the lognormal distribution determines 
the porosity of the medium and is essential to the analysis of 
the transmission of the Lyc. The fit results will be used for the 
modeling of the DIG and discussed in later sections. 

3. MODELS 

3.1. Basic Assumptions 

We employed the "on-the-spot" approximation, in which 
ionizing photons emitted by recombinations to the ground 
state of hydrogen in the ionized gas are assumed to be re- 
absorbed very close to their point of origination. However, 
strictly speaking, the approximation may not be applicable 
to the optically thin case that is the primary concern of the 
present study. The probability of a photon being remitted as 
an ionizing photon is the ratio of the recombination coeffi- 
cient to the ground state of hydrogen (ai) and that to all lev- 
els (a a)- The ratio is 0.3 3-0.43 for a gas in a tempe rature 
range of 5000 - 20, 000 K (iQsterbrock & Ferlandll2006h . This 
portion of the diffuse ionizing photons may escape the sys- 
tem or be reabsorbed elsewhere in the system. Therefore, the 
"on-the-spot" approximation would overestimate the total Ha 
intensity by a factor of < 2 in the optically thin case. The Ha 



morphology would also be probably twice sharper than that 
of a more accurate calculation. These two limitations should 
be kept in mind when interpreting our results. 

The purpose of this study is to simulate the situation in 
which the neutral hydrogen gas has a vertical scale height of 
h, and a density structure that can be described by an absorp- 
tion coefficient hz(z) in the global sense so that some of the 
neutral hydrogen is ionized by the Lyc radiation. Since we as- 
sume that the Lyc photons absorbed in a volume are balanced 
by the immediate recombination, the Ha surface brightness at 
each pixel is directly proportional to the absorbed number of 
Lyc photons within the volume and is independent of the un- 
resolved cloud density structure, such as dumpiness or poros- 
ity. The amounts of Lyc intensities, which are emitted from all 
H II regions and absorbed at a point in the DIG, are co-added 
to calculate the Ha intensity at the point. 

3.2. Thin Slab Models 

We first calculated two-di mensional thin sla b models us- 
ing the method described in Zurita et al.' (2002), not only in 
order to compare our results with theirs, but also to obtain ba- 
sic ideas on the parameter ranges needed to calculate more 
elaborated models. The H I layer in the galactic disk is as- 
sumed to be a thin slab of a constant thickness AH = 2h 
and to have a constant absorption coefficient (kpc~^). We 
took an average H I scale height of h = 100 pc. We cal- 
culated the total radiation field absorbed within a volume 
Ay = (Ar)^A// corresponding to a pixel of the Ha image 
of the face-on galaxy. Here, (Ar)^ is the area correspond- 
ing to a pixel in the image, and the area of the volume sub- 
tended by the source AA = AHAr. Because of geometrical 
dilution and absorption, the transmitted flux of stellar radia- 
tion at distance r from a source is L^e'^ /47Tr^, where Lq is 
the total number of ionizing photons emitted per second at 
the source, and the optical depth r = J ndr. The stellar lu- 
minosity absorbed over volume AV after traveling the optical 
depth r is then given by ^~'^ArLo(AA/47rr^), or, in general, 
^-^(1 -^-^^)Lo(AA/47rr2), where Ar = i^Ar. 

Various absorption coefficients are examined to account for 
the effect of density inhomogeneity. We consider each pixel 
of H II regions as a point source with the effective Lyc lumi- 
nosity equal to that corresponding to its measured Ha inten- 
sity. Given a constant absorption coefficient we obtain the 
Ha surface brightness at pixel /: 

where /scale is the ratio of the total luminosity of DIG to that of 
H II regions. A large fraction of the total Ha emission (DIG + 
H II regions), ~ 0.5, is contributed by DIG, implying that the 
total Lyc luminosity emitted from the sources is about twice 
higher than that absorbed in H II regions. The luminosities 
from the DIG and H II regions are then approximately equal to 
each other (/scale ^1). Here, Vij is the distance from a source 
pixel j to the pixel /, and Tij = nvij is the optical depth along 
the path from the pixel / to / The conversion factors, the Lyc 
to Ha luminosity and the luminosity to flux, are canceled out 

on both sides of the equation. 

It should be noted here that iZurita et al.l (l2002t) ignored the 
factor AA(l - e''^^'). The effect of the factor 1 - e''^^' 
tvAr as K,Ar <C 1) is to lower the amount of produced Ha 
emission as the medium becomes increasingly transparent to 
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the Lyc. The Ha brightness increases with the disk thickness 
Aif . The factor /scale is allowed to vary to take into account 
a larger scale height as well as to match the observed data, 
although it is initially meant to be the ratio of Lyc leakage 
toward DIG to that absorbed in H II regions. In other words, 
the scale height needed to be in accord with the observed data 
ish^ 0.1 /scale kpc. In this way, the scale height of the galaxy 
can be inferred. 

We have varied the effective absorption coefficient from 
K,= 10~^ to 10^ kpc~^. Typical results are shown in Figure 
[2l The scaling factors, defined by /scale = E^data/E^modei, 
were multiplied to the models to match the overall Ha bright- 
ness level. We also show the intensity distributions of the 
model, observed data in the second row, and correlation plots 
between the model and observed flux in the form of two- 
dimensional histograms in the third row. 

We note that, at the optically thin limit of < 0.1 kpc~\ 
as the attenuation factor e~'^^ ^i:^ 1 for a system size of ~ 10 
kpc, the morphologies of the models are basically the same 
as in the case of no absorption (hz = 0) and are determined by 
the geometrical dilution only. In the limit, however, we ob- 
tain much lower Ha intensity than the observed data although 
the overall morphology seems, at first glance, to resemble the 
observed image. The model and observed intensities are, in 
a general sense, linearly proportional to each other, as shown 
for Hi ^ OA kpc~^ in the third row of Figure [2l However, the 
model images are a bit broader near H II regions than the ob- 
served image, as also noted by Zurita et al. ( 2002). As k, in- 
creases, the overall morphology begins to deviate from the ob- 
served one. For large absorption coefficients (e.g., for = 3. 16 
kpc~^), the Ha emissions are mostly concentrated near H II re- 
gions. The correlation between the model and observed data 
also begins to deviate from the linear proportionality, as can 
be noted in third row of Figure [2l We, therefore, fitted the 
correlation between the observed and model intensities with a 

power-law /data = /scale (^model)^- 

The results are shown in Figure O together with /scale- The 
third row of Figure [2] shows the best fit curves in red as well 
as the one-to-one curves in blue. In FigureOb), the slope P is 
a little smaller than 1 even at z^: <C 0.1 kpc~^. It is also shown 
that the slope P decreases rapidly with increasing ksls hz>0.l, 
implying that the observed morphology cannot be reproduced 
well in this range because of fast attenuation of the Lyc while 
it is being propagated into the diffuse ISM. However, when 
the correlation histograms are compared with blue one-to-one 
curves, it is found that the correlation slope is close to 1 up to 
K, = 0.2 kpc~^. Meanwhile, in the limit of < 0.1 kpc~\ the 
scale factor increases with decreasing hz, as shown in Figure 
Oa). Therefore, the optimal solution is obtained with tz ^ 
0.1-0.2 kpc~^. We also note that large scaling factors of ~ 
13-22 are required for these values, implying h ^ 1-2 kpc 
for the generally accepted Lyc leakage fraction of ~ 0.5. The 
required scale height would be ~ 1 - 4 kpc, considering that 
the present results may be overestimated by up to a factor of 
2. These values are too large to be reconciled with the H I 
disk scale height, but consistent with the DIG scale heights 
of ~ 1 - 3 kpc, n ieasured from ou r Galaxy and external edge- 
on ga laxies re.g.. lRevnoldsll989l: Rand 1996; Collins & Randl 
l200lh . 

A serious problem encountered in the present results is an 
extremely low absorption coefficient. Assuming the same H I 
column density of A^hi = 8.6 x 10^^ cm~^ t hrough the galactic 
disk ofh= 100 pc as IZurita et al.l (l2QQ2h . we obtain the ef- 



fective cros s-section (Jeff = /^/( Mii/A//) = 2.3 x 10 ^^(tz/OA 
kpc~^) cm^. IZurita et aH (l2002b obtained the best approxima- 
tion to the morphology of the observed DIG surface bright- 
ness with (Jeff = 5.3 X 10~^^ cm^. This value corresponds 
to tz = 0.2 kpc~\ consistent with our result oftz^ 0.1-0.2 
kpc~^. Assuming a smaller column density of A^hi = 10^^ 
cm~^, we obtain aeff = 2 x 10~^^(hz/0A kpc~^) cm^. These 
cross-sections imply virtually no absorption of the Lyc, i.e., 
as low as ~ 10~^- 10~^ of the photoionization cross-section 
6.4 X 10-1^ cm^ at 912A dVerner & Yak3^[T995h . We also 
note that deff is even lower than the dust extinction cross- 
section of 2.3 X 10-2^ cm^/H at 912 A (Drame 2003). 

3.3. Exponential Disk Models 

For exponential disk models, we employ a three- 

dimensional Mon te Carlo radiation transfer code 

(IWood & Revnoldsl [1999) that was originally developed 
for the simulation of dust- scattering. A modified version has 
been used in the radiative transfer na odels of dust- scattering 
in the Ophiuchus region (|Lee et al.l 12008). The code was 
modified to simulate photoionization in the context of the 
"on-the-spot" aproximation. Lyc photons are emitted from 
point source H II regions, as in thin slab models. We tracked 
the propagation of the Lyc photons from sources as they 
were absorbed and reemitted as Ha photons. The absorbed 
photon was 100% reemitted isotropically as Ha photon in 
the present simulations, while in dust- scattering simulations 
the photon was reemitted with a probability of < a < 1 
toward a new direction determined by a Henyey-Greenstein 
scattering phase function. In other words, our simulations 
were performed by assuming the albedo a = I, asymmetry 
factor g = 0, and a single scattering. The reemitted Ha trav- 
eled toward the observer and was added into the Ha image. 
Absorption coefficients in each density grid were calculated 
assuming an exponential function, k,(z) = tzocxp {-\z\/h), 
where the absorption coefficient hzo (kpc~^) at the galactic 
plane and the vertical scale height h are varied. 

We first confirmed that the Monte Carlo models with /z <C 1 
kpc are consistent with the above thin slab models, and their 
total Ha intensities increase linearly with the scale height. At 
an optical thin limit of /^o <C 1 kpc~\ the morphologies of 
the models are basically the same as the case of no absorp- 
tion (hZQ = 0), just as in thin slab models. Sample models are 
shown in Figure (H together with the intensity distributions 
of the models, observed data in the second row, and correla- 
tion plots in the third row. The scaling factors, defined by 
/scale = E^data/E^modei, wcre multipHed by the models to 
match the overall Ha brightness. Unlike thin slab models, the 
total Ha intensity is saturated at /z ~ 1 / a>:o as the scale height 
increases since the amounts of ionizing photons at high \z\ de- 
clines rapidly and the Ha emission is dominated by the con- 
tributions from low \z\ ^ 1/a^o- Therefore, the scaling factor 
/scale scaturates as h increases. The trend is shown in Figure 
[5] for the models with hcq = 0.4, and 10.8 kpc~^ As hcq in- 
creases, the Ha morphology becomes strongly concentrated 
near the H II regions, and the intensity distribution becomes 
wider as its low-intensity side is elongated. The elongation 
of the low-intensity side of the intensity histogram is shown 
in the blue curve of Figure Hfd), in which no /scale was mul- 
tiplied. For a given absorption coefficient, the mophological 
coincidence of the thick disk models with the observation is 
better than that of the thin slab models, as the Ha emission 
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from the thick disk models are from more extended regions. 
In other words, the intensity distribution is narrower than that 
for the thin slab models. This can be easily noticed by com- 
paring the model of hZo = l kpc~^ in Figure[2land that of hZo = 2 
kpc~^ in Figure (H Even with a higher absorption cofficient, 
the thick disk model is less concentrated near the H II regions 
than the thin slab model. As for thin slab models, we fitted the 
correlation between the observed and model intensities with 
a power-law and plotted the best fit curves in red lines and 
one-to-one curves in blue lines in the third row of Figure HI 
We found that the models with hZQ = 0.4-0.8 kpc~^ predict the 
observed data nicely, and the scale height required to explain 
total Ha intensity is ~ 1 - 2 kpc~^ . 

The problem of the extremely low absorption rate remains 
unsolved with the exponential disk models, although four 
times higher hzo are obtained. Assuming the same H I column 
density over 0.1 kpc, as for the thin slab models, we obtained 
the effective absorption cross-section (Jeff = (0.9-1.9) x 10~^^ 
cm^, as low as ~ 10~^ of the photoionization cross-section. 
The obtained scale height is consistent with that from thin slab 
models. We note again that the scale height would be ~ 1-4 
kpc, considering that the present results may be overestimated 
by up to a factor of 2. 

4. OPACITY OF THE TURBULENT ISM TO THE LYC 

We now examine whether such a small attenuation coeffi- 
cient obtained in §3 is in accord with the current knowledge 
on the ISM density structure. What we measure from ob- 
servation is the average column density (N) or optical depth 
(r) because of the finite spatial resolution. In the turbulent 
medium showing the lognomal nature of column density, the 
effective optical depth that the radiation field suffers, given an 
average optical depth (r), can be calculated by 

Teff = -ln J e~^P(lnr)dlnr^ 

where the lognormal distribution is 
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P(lnr) = 



1 



27rcrinr 



exp 



(Inr-lnro)^ 



(iFischera. Dopita. & Sutherland! l2003l) . Here, (Inr) = Inro 
and ainr are the mean value and the standard deviation in 
logarithmic space. The mean values in linear and logarith- 
mic scales are related by Inro = ln(r) where (r) 
is the mean in the linear scale. The relative standard devi- 
ation (contrast) in the linear scale, ar/ (r) = crr/{r)^ is re- 
lated to the standard deviation in the logarithmic scale by 
^r/ir) ~ ^^P(^inr)~ ^- relative deviation of optical depth 
crr/{r) is always smaller than the density contrast crp/{p), as 
the accumulation of the density along a path length would 
average out the density contrast ( Fischera & Dopita 2004). 
Given the observational or theoretical value of the column 
density contrast cr^/(r), we can, then, obtain the lognormal 
distribution that describes the column density, and estimate 
the effective optical depth, Teff, of the Lyc using these rela- 
tionships. 

Observationally, the turbulent ISM has been usually 
investigated by a power-law fit to the power spectra of the 
observable quantities proportional to column densities. The 
power-law indexes have been observed tohcy ^ -2.4 ~ -3.7 
from the warm and cold ISMs in our Galaxy and the Small 
Magellanic Cloud (Armstrong. Cordes. & Rickett .198ll: 



2007). The power-law indexes can be related to the density 
contrast cr^/^p^ by combining two relationships between 
the Mach number M and 7, and between M and (Tpiip) - 
We combine the numerical results of Kim & Rvul (12 0051) . 
iKritsuk. Norman. & PadoanI (|2006), and Padoan et al. ( 2o3, 
and find a relation between the Mach number and the 
power-law index of the density power spectrum, 

7 = (-3.81±0.07)M(-0-i™\ 

as shown in Fig ure [6l From this, we estimate M ^ \ - 12 for 
the observed 7. iFischera. Dopita. & Sutherland! (|2003|) found 
that the empirical reddeni ng law for starburst galaxies is con- 
sistent with M^\3-22. iPadoan. Jones. & Nordlundl (Il997h 
found that cFp/{p) ^ 0.5M through numerical simulations. We 
then obtain the upper limit on the column density contrast of 
^t/{t) ^ ^p/{p) ^0.5-11. In Figure [71 we show Teff as a 
function of cFt/{t) for various (r). 

Although, in real situations, the H I column density is an 
order of ~ 10^^- 10^^ cm~^ over 0.1 kpc scale, we calculate 
the effective optical depths for various column densities to 
demonstrate how small an absorption coefficient is required 
by the "standard" photoionization scenario. Assuming H I 
column densities of A^hi = 10^^, 10^^, and 10^^ cm~^ over 0.1 
kpc, we have (r) = A^HiCTphoto (1 kpc/0.1 kpc) - 64,000, 6400, 
and 640 over a distance of 1 kpc, respectively. We then obtain 
effective absorption coefficients of ^ 11 - 159, 6.2-92, 
and 3.5-45 kpc~^ from Figure[7l for cr^/(r) ~ 0.5 - 1 1. Note 
that these values are the lower limits for various H I column 
densities, but are still much larger, by several orders of mag- 
nitude, than the /^o ^ 0.4-0.8 kpc~^ required to account for 
the DIG by the "standard" photoionization mode l. 

Recently, IMiville-Deschenes & Martini (120071) found that 



V(P> 



0.8 for the warm neutral medium based on the kine- 



matic H I observation. 'Hill et al." (2008) showed for the first 
time that the Ha intensity PDF of our Galaxy is well repre- 
sented by a lognormal distribution and found that M = 1.3 - 
2.4 through comparison with the turbulence simulations. We 
performed the same fit, as shown in Figure [TJb) and obtained 
a width of the lognormal distribution, ain/ ~ 0.36 ± 0.03 dex, 
~ 1.5 times wider than that of our Galaxy. We obtain M < 
6.5-7.3 for the DIG of M 51 by comparing the result with 
Table 1 in Hill et al. (.2008 ). Assuming M = 8 and cr^/(^) ^ 4, 
we have much stronger constraints on the absorption coeffi- 
cient of /^eff > 19.3, 10.8, and 6.1 kpc'^ for A^hi = 10^^ 10^^ 
and 10^^ cm~^ over 0.1 kpc, respectively. 

The models calculated with these absorption coefficients 
yield Ha morphologies significantly different from the ob- 
served one. This implies that the Lyc photons leaked from 
the H II regions are immediately absorbed by the diffuse ISM, 
even if it has escaped out of the H II regions. For compar- 
ison, we calculated the Mach numbers required to have the 
effective absorption coefficient /t: = 0.8 kpc~\ assuming that 
the same relation between the standard deviation of column 
density and the Mach numer is still valid in this extreme case: 
M > 2.0 X 10^ 1.8 X 10^, and 1.6 x 10^ for the media having 
A^Hi = 10^^, 10^^, and 10^^ cm~^ over 0.1 kpc, respectively. 

5. DISCUSSION 
5.1. Vertical Profile 
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By comparing the Ha emission data with models, we found 
that the vertical scale height of the face-on galaxy M 5 1 must 
be ~ 1-2 kpc. Even though this value is obtained from 
the problematic absorption coefficient, we note that the scale 
height is consistent with the results not only from our Galaxy, 
but also from exter nal edge-on galaxies (Reynolds 1989; Ran J 
119961: ICoUins & Ra nd 2001). It might also be noteworthy to 
examine how the galaxy M 5 1 may look when the galaxy is 
viewed edge-on. In FigurelH we show the Ha map, assuming 
hZo = OA kpc~^ and h=l.5 kpc, which are expected when the 
galaxy is viewed edge-on, and the vertical profile of Ha im- 
age. The vertical profile of the model with the a^o = 0.4 kpc~^ 
and h=l.O kpc are shown as well for comparision. The Ha 
intensity from H II regions are not shown in the figure. 

It is surprising that the simulated vertical Ha profile can 
be represented by a function consisting of two or more- 
exponential components. Observationally, the profile might 
be recognized by two-exponential components because of low 
surface brightness at high \z\. If the Lyc photons travel only 
perpendicularly to the galactic disk, the Ha profile would 
have the same scale-height as that of the exponential disk. 
However, the Ha radiation observed at a point with verti- 
cal height z consists of the contributions from the Lyc pho- 
tons that have traveled longer distances, i.e., (r^ > \z\, 
where r is the projected distance onto the galactic plane 
from the Lyc source. Therefore, the Ha intensity declines 
faster than that naively expected from the assumed vertical 
height of the exponential disk. The smaller the vertical dis- 
tance, the more important the effect becomes. Beyond the 
height correponding to the system size ~ 10 kpc, the ver- 
tical profile converges to that of the disk, as the distances 
traveled by the Lyc become close to the vertical height, i.e., 
(r^+z^)^/^ ^ \z\. The trend is noticeable in Figure [8tb). We 
then have a vertical profile that resembles two-exponential 
components, within a few kpc. The profiles in Figure [Ub) 
are fitted with 0.19^""/^-^ +0.33^-^/^-^^ for h=l.5 kpc, and 
0.23^-^/0-^0 + 0.32^-^/0-26 for h=l.O kpc, and the fit profiles 
are denoted by diamonds. We have masked the regions at 
|z| < 0.3 kpc in fitting the vertical profile to minimize the ef- 
fect due to the assumption that the H II regions are located in 
the galactic plane (z = 0). The vertical height of the smaller 
exponential component is expected to depend not only on the 
vertical profile of H I gas, but also of the details of the vertical 
extents of the H II regions. 

The Ha vertical emission profiles of edge-on galax- 
ies are known to be well represented by two exponen - 
tial components ( "Rand' '1997'; 'Miller & Veilleux' '2003a'). 
IWang. Heckman. & Lehnert (1997) identified the two compo- 
nents kinematically through the observation of face-on galax- 
ies, including M 51, a "quiescent" DIG with a high ioniza- 
tion state and small scale height (few hendred pc), and a "dis- 
turbed" DIG with a high ionization state and moderate scale 
height (0.5-1 kpc). The origin and evolution of the extended, 
or extraplanar, DIG seen in edge-on galaxies are not clear, but 
two alterna tives are generally considered (e.g., i Heald et aH 
[20Q6ll2007h . The gas is generally considered to participate in 
a star formation-driven disk-halo flow, su ch as that described 
in the fountain or chmi ney models ( Sha piro & Fi eld 1976); 
iNorma n & Ikeuchi"1989). Alternatively, the gas could have 
been accreted from t he intergalactic medium or fro m compan- 
ion galaxies (e.g., Ivan der Hulst & Sancisill2005h . The ver- 
tical structure seen in our models seems to support the star 
formation-related origin of the extraplanar DIG. However, 



without kine matical information, as was used in iHeald et aH 
(120061 l2007h . it is not easy to examine the two-components 
predicted in these models with = 0.4 kpc-^ do correspond to 
the kinematically identified components, and/or to determine 
how much of the DIG is associated with the quiescent com- 
ponent and how much is contained in extended disturbed lay- 
ers. In fact, it is not only surprizing, but also very puzzling, 
that the models with a problematic absorption rate predicts 
not only the apparent morpologies of M 5 1 , but also the two 
exponential components seen in edge-on galaxies. 

5.2. Density Stucture of the ISM 

As suggested by numerical s imulat ions (e.g., 
IVazquez-Semadenil 119941: iKim & Rvul I2005h . the den- 
sity distribution of turbulently disturbed ISM is represented 
by a lognormal distribution, of which the width depends on 
how the medium is turbulent. Padoan, Jones, & Nordlundj 
(1997) found that the observational result of ILada et al] 
(1994), which was that the variation of the stellar extinction 
in dark clouds correlates with the mean extinction, is con- 
sistent with the logno rmal distribution of the dust density. 
IWada. Spaans. & Kiml ([2000) showed that the luminosity 
function of the H I column density in the Large Magellanic 
Cloud is lognormal. More recently. Hill et al. (2008) found 
that the Ha emission measure perpendicular to the Galactic 
plane from the DIG in the Milky Way is well fitted by a 
logn ormal distribution, with a stand ard deviation of crin/ ~ 0.2 
dex. Berkhuijse n & Fletched ([2008) showed that the PDFs of 
the average density, not only of the DIG, but also of the H I 
gas in our Galaxy, are close to lognormal, with a standard 
deviation (T(„) ~ 0.3 dex. In Figure (Hb), we showed that 
the PDFs of the Ha intensity, i.e., emission measures, from 
the DIG and H II regions in M 5 1 are well-represented by a 
lognormal function. The standard deviation obtained from 
the DIG in M 5 1 is a bit higher than those of our Galaxy. 

Wang , Heckma n, & Lehnert (1997) studied kniematics of 
the DIG in M 51, together with other galaxies and found 
that the DIG is more disturbed kinematically than the gas 
in the giant H II regions. We, thus, compare the widths of 
Ha intensity distributions of the DIG and H II regions in Fig- 
ure [lib) and found no significant difference between them, 
except an extended tail toward low intensities shown in the 
DIG only. However, it is unclear whether the extended tail 
might be related to the more distu rbted component seen by 
IWang. Hec kman. & LehnertI (Il997h . 

In §4, we found that an unrealistically small absorption, 
even smaller than dust-extinction, is necessary to explain the 
global features of the DIG with the Lyc leakage from the H II 
regions. The turbulent ISM with lognormal density PDF can 
be much more transparent to the Lyc, up to an order of 1 ~ 3 
depending on the width of the PDF, compared to homoge- 
neous ISM, as shown in Figure [71 However, such a small 
attenuation found here cannot be reconciled with these val- 
ues, as discussed in §4. In addition, the attenuation rate is 
not consistent even with the dust optical depth of Tdust ~ 13 at 
912A (Ay ~ 3) obtained by comparing Ha a nd Paa emissions 
measured to ward the H II regions of M 5 1 ( Quillen & Yukital 
120011: IScovil le et al. 2001). Such large amounts of dust sur- 
rounding H II regions should absorb Lyc photons immediately 
in the vicinities of the H II regions. 

Radiative transfer models have been successful in explain- 
ing the dust scattering and absorption for our Galaxy and 
external galaxies with reasonable amounts of dust, consis- 
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tent with the observed amounts of int erstellar dust (e.g., 
ISchiminovich et al.ll2QQTl: lBianchill2007h . Unless the topol- 
ogy of the interstellar dust and neutral hydrogen are com- 
pletely different and decoupled from each other, the radia- 
tive transfer of the Lyc should not differ significantly from 
that expected from observed amounts of neut ral hydrogen. In 
their study of dust scattering in clumpy media. Iwitt & GordonI 
([1996 ) showed that the effective optical depth of a system 
is mainly determined by the covering factor, i.e., the frac- 
tion of 47r steradian of a solid angle filled by clumps as 
seen by the source rather than by the individual optical depth 
of clumps. The transmission of the Lyc is also mainly de- 
termined by the covering factor of dilute regions seen by 
the source. Such a large covering factor of low density re- 
gions may be achieved by fountain flow or chimne ys around 
OB associations ( Shapi ro & Fieldl 119761: iNorman & Ikeuchi 
[1989). Filaments seen in edge-on galaxies may originate from 
the star formation-driven structures. Meanwhile, Rossa et al. 
(2004) presented high spatial resolution (5 pc) observations 
of the extraplanar DIG in the edge-on galaxy NGC 891 and 
found no clear broadening of the filamentary structures to- 
ward high \z\ expected from individual chimneys. They, there- 
fore, concluded that the chimney scenario is most likely not 
responsible for the gas and energy transport into the halo. 
In addition, iWang. Heckman. & Lehnert (1997) found that 
M 51 has a pronounced "diffuse" Yla emission based on 
its high value of m ean/rms of the Ha surface brightness. 
iDove. ShuU. & Fer rara ( 2000) showed that the dense swept up 
shell of material by the dynamical action of supernovae and 
stellar winds would trap ionization photons within superbub- 
bles. Given the diffuseness or smoothness of the Ha emission 
and the approximate equality of Ha luminosities of the DIG 
and H II regions, the ISM should have an extremely unreal- 
istic configuration, completely evacuated straight thin holes 
extending over the order of 1 kpc, randomly covering ~ 1/2 
of the 47r solid angle seen by the sources. 

The atomic ISM may be divided into two phases, except 
the hot gas of temperature ~ 10^ K: the neutral hydrogen 
gas (or clouds) and the DIG (or interclouds) phases. Our 
aim was to investigate the situation in which the stellar Lyc 
radiation propagates into the "initially" neutral ISM, ionize 
some portion of the gas, and maintain two phases. Since the 
ISM is not fully ionized, the H I gas prevents the propaga- 
tion of the Lyc photons to maintain the ionization state of the 
DIG. However, it would be worthwhile to imagine an ISM in 
which the H I gas phase has no effect on the propagation of 
the stellar Lyc photons and the freely transmitted Lyc radia- 
tion maintains the ionization state of the DIG. In other words, 
assume that the H I gas phase is completely evaculated some- 
how, and the DIG is left alone. Then, the radiative tranfer 
in the medium will be determined by the observed property 
of the DIG that is deduced from the Ha intensity distribu- 
tion in §4. Adopting M = 8, crr/{r) = 4 from the DIG of M 
51, and A^hi ^ 1 x 10^^ cm^ ov er ^ 1 kpc estimated in the 
Galactic plane (lReynoldslll991h . the effective absorption co- 
efficient is ~ 10.8 kpc~^. The model for these parameters is 
shown in Figure Hfd). The model shows a highly concen- 
trated Ha profile near the H II regions and predicts only about 
1 //scale ~ 26% of the observed Ha luminosity. Since the rel- 
ative deviation of the column density is an upper limit, the 
amount of Ha that can be produced in this medium is even 
smaller, i.e., < 26%. A more stringent constraint would be 
obtained, assuming a Mach number of M =^ 1.2 -2 A esti- 



mated from our Galaxy (iHiU et al.ll2QQ8l) . Therefore, the OB 
associations in the H II regions are not likely to be a main 
source in ionizing the DIG. The only way to resolve this prob- 
lem would be to assume alternative ionizing sources outside 
of H II regions, in-situ of the DIG. If this is the case, the sim- 
ilarity in the Ha morphologies between the models and the 
observational data may be a coincidence. 

6. CONCLUDING REMARKS 

We investigated whether the global features of the DIG can 
be explained in the context of the "standard" photoionization 
model. The most important result of this study is that we 
quantified the topology or "porosity" of the ISM required by 
the "standard" photoionization scenario in the contexts of the 
lognormal PDF. We found that even without the H I gas phase, 
the DIG phase cannot be maintained by the stellar Lyc radia- 
tion alone, implying that main ionizing source is not the OB 
associations in the H II regions. Since the Lyc that leaked out 
of the H II regions can account for, at most, only < 26% of the 
Ha luminosity of the DIG, the alternative sources should play 
a significant role in ionizing the DIG, especially in the inter- 
arm regions. Even with the mophological coincidence of the 
Lyc leakage models with the Ha data, we cannot rule out the 
possibility that alternative ionizing sources are generally cor- 
related with star formations and, thus, with the H II regions, as 
various gas and dust components are well correlated with each 
other, and the sources are playing a significant role in ioniz- 
ing the DIG. The strong correlation between the DIG and H II 
regions does not necessarily indicate the direct association of 
the DIG with the H II regions. 

Other authors have arrived at similar conclusions, both 
in terms of producing the amount of ionized gas observed 
and in producing the observed spectral signature of the 
DIG, that alternative sources other than hot stars in the 
H II regions are needed in heati ng and/or ionizing the DIG. 
[Rand. Kulkarni & Hestel (119901) estimated that the total en- 
ergy input by chimneys into the halo to account for the DIG 
of NGC 891, a galaxy that in many ways is very similar to 
our own, is a factor 10-20 higher than the estimates for our 
Galaxy. Although the supernova rate is likely to be a few 
times higher than in our Galaxy, the required superbubble for- 
mation rate was inconsistent with the observations. 

The emission-line ratios of the DIG are significantly 
different from those of typical H II regions. In the DIG, 
[Si II] A6716/Ha and [N II] A6583/Ha are enhanced 
relative t o the H II re gions (Haffner, Revn olds. & Tufte 
T999I: iFerguso n. Wvse , & Gallagher 199^ iRand 
T997I: Greenawalt ltair Il998l: ICo llins & Rand 2001; 
Hoopes & WalterbosI 120031: IVoges & Walterbos. i2006) . 
The pure photoionization models in diluted media predict 
the enhancement very well. However, the models predict 
that the ratio [Si II]/[N II] should increase with z, in con- 
trast to the observations that the ratio versus the distance 
from the H II regions rema i ns ne a rly constant in the DIG 
(iHaffner. Revnolds. & Tuft3 [T999I: ICollins & Randl I200TI) . 
The behaviors of these low ionization lines can be explained 
by the gas temperature increase alone, without having a 
secondary source of ionization (Haffner. Revnolds. & Tuftd 
[T999h . On the other hand, the ratio [O'lII] A5007/Ha 
cannot be fully explain ed by a temperature increase alone 
(IColhns & Randl 120011) . According the "standard" pho- 
toionization scenario, the ratio [O III]/Ha should decreases 
with z or the distance from the H II regions in a pho- 
toionization scenario. However, the opposite trend or the 
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relative cons t ant with z is seen not o nly in edge-on galaxies 
(iRandl 119981: ICoUins & Randl l2001h. but also in face-on 
galax ies including M 51 (IWang. He ckman, & Lehnei 
19971: iHoopes & WalterbosI I2Q03E Fvoges & Walterbo 
20061) . One possible way to explain this high line ratio 
is to assume a secondary ionization source that con- 
tributes a small but in creasing fraction of the Ha with z. 
ICollins & Randl (1200 11) explored the p ossibility that ioniza- 
tion by shocks j^ ull & McKee Jl979h o r turbulent mixing 



layers (Slavin, Shull, & BegelmanI Il993h is responsible for 
some of the DIG emission of the edge-on galaxies NGC 
4302 and UGC 10288. They found that the composite 
models, of which up to ~ 16-45% of the Ha arises from 
shocks, better reproduce the variations of line ratios than the 
photoionization models alone, although no model reproduces 
the line ratios to within the errors. Similar conc l usions 
were obtained for other edge -on galaxies by Rand" ("1998) 
and Miller & Veilleuxl (!2003bl). ^Hoopes & Walterbos (2003) 
tested the feasibility of "leaky H II region" models, in wich 
the transmitted ionizing continua of density-bounded H II 
regions power the DIG, and concluded that other ionization 
or heating mechanisms play a role in ionizing the DIG. In 
addition, the He I A5876/Ha line intensity ratio in the DIG of 



our Galaxy is found to be significantly less than that observed 
in H II regions, challen ging the O star ionization models 
(iRevnolds & Tuftd[l995h . 

Ano t her scenario investigated by iHoopes & WalterbosI 
(|2000|),lHoopes, Walterbos, &£othun| (12001) is that field OB 
stars ionize some of the DIG. They measured the Ha/FUV 
intensity ratios of ten spiral galaxies, and concluded that late 
OB stars in the field are indeed impo rtant ionization sources. 
However. I Voges & WalterbosI (|2006|) concluded that field star 
ionization models do not fit the data of M 33 well and that 
about 30% leakage from H II regions is necessary. Accretion 
from the intragalactic medium or companion galaxy would 
be an att ractive alternative for the origin of, at least some of, 
the DIG. iHeald et al.l (I2006L 120071) compared the kinematics 
of optical emission lines obtained from two edge-on galaxies 
NGC 891 and N GC 4302 with an balhstic m odel of a galactic 
fountain model (Collins, Benjamin, & Rand"2002) and found 
a strong disagreement between the observed kinematics and 
those predicted by the model. As they noted, galaxies may 
reside in halos populated by a various mixture of accreting 
gas and star formation-driven gas. However, further analyses 
are required to quantitatively estimate their contributions in 
ionizing the DIG. 
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Fig. 1 . — (a) ¥La image of M 51 in logarithmic scale (EM, cm ^ pc). The H II regions were masked, (b) Histograms of Ho; intensities for the DIG and H II 
regions. Blue curves denote the best-fit lognormal functions, as discussed in §4 and §5. 
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Fig. 2. — Thin slab models calculated for various absorption coefficients. Predicted model images with the same color scales as in Fig.[T] intensity histograms 
of the model and observed data, and correlation plots between the model and observed data are shown in the first, second, and third rows, respectively. 
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Fig. 3.— Best fit parameters with /^ata = /s*aie (Anodei)^- The dotted fine was obtained by /scale = Z)4ata/ Z)Anodei- 
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Fig. 4. — Exponential disk models calculated for various absorption coefficients and scale height: (a) k,o = 0.1 kpc~^ and h = 3 kpc, (b) kq = 0.4 kpc~^ and 
h = l.5 kpc, (c) ^co = 2.0 kpc~^ and h = 1 kpc, and (d) kq = 10.8 kpc~^ and h = l kpc. Predicted model images with the same color scales as in Fig. [T] intensity 
histograms of the model and observed data, and correlation plots between the model and observed data are shown in the first, second, and third rows, respectively. 
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Fig. 5. — Variation of the scale factor, /scale = X^^data/ X^^modeh for the exponential disk model with the absorption coefficient = 0.4 kpc ^ and 10.8 kpc ^ . 
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Fig. 6. — Power-law slope versus the Mach number from numerical simulations of the turbulent ISM, and the best-fit curve. 
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Fig. 8. — (a) Model for the galaxy M 51 when it is viewed edge-on, assuming kq = 0.4 kpc~^ and h = 1.5 kpc, (b) the vertical profile of the Ha intensity for 
two models: kq = 0.4 kpc~\ h=1.5 kpc, kq = 0.4 kpc~\ and h=1.0 kpc. Here, the diamond symbols denote two-exponential fits for |z| > 0.3 kpc. 



